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Abstract 

We have employed the constrained variational method to study the influence of spin polarization 
on the ground state properties of liquid ^He. The spin polarized phase, we have found, has stronger 
correlation with respect to the unpolarized phase. It is shown that the internal energy of liquid 
^He increases by increasing polarization with no crossing point between polarized and unpolarized 
energy curves over the liquid density range. The obtained internal energy curves show a bound 
state, even in the case of fully spin polarized matter. We have also investigated the validity of 
using a parabolic formula for calculating the energy of spin polarized liquid ^He. Finally, we have 
compared our results with other calculations. 
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I. INTRODUCTION 



The spin polarized liquid ^He is an interesting quantum many-body system which can be 
experimentally examined. In fact, it is expected that this phase of liquid ^He has a large 



life time to be observed in a quasi-thermodynamic equilibrium 




a. 



Some theoretical investigations have been done for spin polarized liquid ^He using differ- 
ent approaches such as Green's function Monte Carlo (GFMC), Fermi hyper-netted chain 



fly , 3, Q, Q, Q, Q. 



(FHNC), correlated basis functions (CBF) and transport theory 
Recently, we have studied the unpolarized liquid ^He and calculated some of its thermo- 
dynamic properties at finite temperature. A good agreement between our results and cor- 



responding empirical values has been shown 



14 



151]. In these calculations, constrained 



variational method based on the cluster expansion of the energy functional has been used. 



This method is a powerfu 

matter Q, 0, Q, Q, y , 



microscopic technique used in many-body calculations of dense 



21 



22l |. In this article, we employ this method to investigate the 
ground state of spin polarized liquid '^He. 



II. CONSTRAINED VARIATIONAL CALCULATION OF SPIN POLARIZED 
LIQUID ^He 

To calculate the ground state energy of spin polarized liquid '^He, consisting of in- 
teracting atoms with A^'-"'"^ spin up and N^~^ spin down atoms, = A^'-"'"-' + A^^"-*, we use 
the variational calculation based on the cluster expansion of the energy functional 2J. We 
consider up to the two-body energy term in the cluster expansion, 

E = Ei + E2. (1) 

The one-body energy per particle for the spin polarized matter is given by 



E, = + e[-^ 
3 /n'^ 



10\2mJ^ ' 



(2) 



^;i+^)5/3 + (i_^)5/3 

where p is the total number density, 

p = pW+p(-), (3) 



2 



and the spin asymmetry parameter, ^, is defined to be 

iV(+) - iV(-) 
? = - 



N 



(4) 



^ shows the spin ordering of the matter which can get a value in the range of ^ = 0.0 
(unpolarized matter) to = 1.0 (fully polarized matter). 

To obtain the two-body energy term, E2, we can start from the known equation 

1 



^2 = ^E<^jW(12)Kj>a, 



where 



W{12) 



m 



V/(12) +f\l2)V{12), 



(5) 



(6) 



and V{12) is the two-body potential between the helium atoms (we use the Lennard- Jones 
potential with e = 10.22 K and a = 2.556 ^4) and /(12) is two-body correlation function. 
By considering the |z > as a plane wave and performing some algebra, we have derived the 
following relation for the spin polarized matter. 



1 + 0'^'(4^Vi2) + (1 - 0"<i'\kPn2 



where 



(X 



(^)^3 



Wiru), (7) 



(8) 



sin(2;*^*^) — x^*^ cos(a:^*'') 

and x*^*^ is kp'^ru or kp'^ru- kp'^ = {Qn'^p^^'>y^^ and kp'^ = {Gn'^p^'^y^^ are the Fermi 
momentum of spin up and spin down states, respectively. By minimizing the Eq. ([7j) with 
respect to the /(12), the following differential equation is obtained: 



d_ 
dr 



L(r)/'(r) 



m 



V{r) + X 



where 



L(r) 



i + 0'^'(4^'ri2) + (i-0'^'(^F^ri2 



L(r)/(r) = 0, 



(9) 



(10) 



and A is the lagrange multiplier which imposes the constraint of the two-body wave function 
normalization. Eq. iQ can be solved numerically and the correlation function /(12) and 
finally the internal energy of system are obtained. 
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III. RESULTS 



The numerical results for the two-body correlation function of liquid ^He are given in 
Fig. [TJ This shows that for the higher value of spin asymmetry parameter, the correlation 
function increases more rapidly and reaches to the limiting value (/(r) = 1) at smaller value 
of r. From Fig. [H we can see that in the case of higher polarization, the helium atoms have 
a stronger correlation at short relative distance. 

The numerical values of kinetic, potential and internal energy of liquid ^He for the different 
polarizations versus total number density are presented in Tables Hi ll VI Our results of internal 
energy for different spin asymmetry parameters are compared in Fig. [21 This figure indicates 
that as polarization increases, the internal energy gets the higher values over the liquid 
density range and there is no crossing point between the internal energy curves of polarized 



and unpolarized cases in this region. This behavior is in agreement with experiment 10|, 112, 



13|. A comparison between Tables Ht [IT | IIIII and IIVI shows that for each density, the kinetic 



energy (potential energy) increases (decreases) by increasing spin asymmetry parameter. 
Over the liquid density range, the increasing of kinetic energy dominates which leads to 
the higher internal energies. From Fig. 2, it can be also seen that for all values of the 
energy curve has a minimum which shows the existence of a bound state for this system, 
even for the fully polarized matter = 1.0). We see that by increasing the spin asymmetry 
parameter, this minimum point of energy curve shifts to the higher densities. 

To compare our method with the well-known many-body techniques, we present the 
results of Green's function Monte Carlo (GFMC) and Fermi hyper-netted chain (FHNC) 
calculations in Figs. [3B In Fig. [31 the internal energy of fully polarized and unpolarized 
liquid ^He calculated with GFMC method 12|, |l3| have been compared with our results. 



There are two points that we can mention from this figure. First, the crossing point prob- 
lem exist and secondly, for all densities approximately grater than 0.014A~^, the energy of 
polarized case is lower than unpolarized case which are not acceptable from the experiment. 
However, these problems do not exist in our method, although we should include the three- 
body correlation effect to obtain a better result. In Figs. [1}[H1 the results of FHNC method 
for different choices of wave function 10|| are presented. We can see from these figures 



that before considering the backflow effect (momentum dependent two-body correlation), 
the polarized curve is always lower than the unpolarized curve. However after performing 
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extra calculations, the appropriate results have been obtained in Fig. [61 This shows that 



our constrained variational method can do the job much sirap 



22 



er, although we still need to 



24|. 



add the three-body cluster energy to obtain a better result 

There is a similarity between spin asymmetry parameter in our calculations and isospin 
asymmetry parameter in nuclear matter calculations in which the energy of asymmetrical 
nuclear matter can be calculated using the parabolic approximation. In this approximation, 
one considers only the quadratic term in asymmetry parameter as well as the energy of 
symmetric matter 



2l| . In a similar way, we can define the following relation for the internal 
energy of spin polarized liquid ^He: 

E{p, = E{p, e = 0.0) + aasymip) i\ (11) 

where the above equation gives the definition of spin asymmetry energy aasym. as: 

aasymip) = E{p, e = 1.0) - E{p, e = 0.0). (12) 

£'(p, ^ = 0.0) is the internal energy of unpolarized matter which has a symmetric configu- 
ration in spin state and E{p, ^ = 1.0) is the internal energy of fully polarized matter. It is 
now interesting to see if this approximation, Eq. fllip . agrees with our microscopic calcula- 
tions. For this purpose, we have compared the internal energy of spin polarized liquid ^He at 
different ^ for both microscopic calculation and using parabolic approximation in Tables IVl 
andlVH These tables indicate an agreement between these two approaches, specially at high 
densities. 



IV. SUMMARY AND CONCLUSION 

We have considered a system consisting of Helium atoms (^He) with an asymmetrical 
spin configuration and derived the two-body term in the the cluster expansion of the energy 
functional. Then, we have minimized the two-body energy term under the normalization 
constraint and obtained the differential equation. The numerical results of internal energy 
of this system have been presented for different values of spin asymmetry parameter and 
density. It is found that as the polarization of liquid ^He increases, the two-body correlation 
becomes stronger. Over the liquid density range, our results show that the internal energy 
of liquid ^He increases by increasing spin asymmetry parameter with no crossing point be- 
tween polarized and unpolarized energy curves. This shows an agreement with experimental 

5 



results. It is also seen that there is a bound state for all values of polarization. The validity 
of using parabolic approximation in calculating the energy of spin polarized matter is shown. 
Therefore, by using this approximation one can preform calculations for the spin polarized 
matter much simpler, we have also compared our results with other calculations to show 
that why of our constrained variational method is a powerful technique. 
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TABLE I: Kinetic energy (KE), potential energy (PE) and internal energy of unpolarized liquid 
^He = 0.0) versus total number density. 



Density fA~"') 


KE (K) 


PE (K) 


Internal Energy (K) 


0.003 


0.960 


-0.832 


0.128 


0.005 


1.349 


-1.578 


-0.228 


0.007 


1.689 


-2.341 


-0.651 


0.009 


1.997 


-3.176 


-1.178 


0.011 


2.283 


-3.864 


-1.581 


0.013 


2.552 


-4.049 


-1.496 


0.015 


2.807 


-3.386 


-0.578 


0.017 


3.052 


-1.422 


1.629 


0.019 


3.287 


2.321 


5.608 



TABLE II: As Table H but for spin polarized liquid ^He at ^ = 1/3. 



Density (A~^) 


KE (K) 


PE (K) 


Internal energy (K) 


0.003 


1.019 


-0.922 


0.096 


0.005 


1.433 


-1.701 


-0.268 


0.007 


1.794 


-2.482 


-0.688 


0.009 


2.121 


-3.256 


-1.135 


0.011 


2.424 


-3.957 


-1.532 


0.013 


2.710 


-4.154 


-1.444 


0.015 


2.981 


-3.505 


-0.523 


0.017 


3.241 


-1.555 


1.685 


0.019 


3.491 


2.176 


5.667 
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TABLE III: As Table H but for spin polarized liquid ^He at ^ = 2/3. 



Dpnsitv ( A 


KE (K) 


PR (K) 


Tnf PTTl 1 PTl PVCAT" / T\ 1 
-LlltCllldl Clld \ -LV f 


yj.Kjyjo 


1 909 


1 1 "^"^ 


OfiQ 


u.uuu 




1 QfiS 
-i.yuo 


97Q 

-U. Zr ( U 






Z. I uu 




OOQ 


9 AQQ 


1^1 R 


1 01 S 


01 1 




A 94n 


1 ^^^9 

-1 .OOZr 


0.013 


3.194 


-4.477 


-1.283 


0.015 


3.514 


-3.863 


-0.349 


0.017 


3.819 


-1.945 


1.874 


0.019 


4.114 


1.754 


5.868 



TABLE IV: As Table H but for fully polarized liquid ^He = 1.0). 



Density (A"^) 


KE (K) 


PE (K) 


Internal energy (K) 


0.003 


1.524 


-1.403 


0.120 


0.005 


2.142 


-2.227 


-0.085 


0.007 


2.681 


-2.967 


-0.286 


0.009 


3.170 


-3.648 


-0.478 


0.011 


3.624 


-4.268 


-0.644 


0.013 


4.051 


-4.850 


-0.799 


0.015 


4.457 


-4.478 


-0.021 


0.017 


4.845 


-2.615 


2.229 


0.019 


5.218 


1.033 


6.251 
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TABLE V: Internal energy (K) of spin polarized liquid '^He versus total number density at ^ 
computed with both microscopic calculation and using parabolic approximation. 



Densitv ('A"''! 


A/Tirrnsronir C^almlatinn 

± V A± v>i. v_/0 V-j\-J L/± v> V_.- (_h± LXJ. CAj Lj ±\-J LL 


ParaHolir A DDrnxiTn a tion 


0.003 


0.096 


0.127 


0.005 


-0.268 


-0.212 


007 


-0 688 

\J .WOO 


-0 61 1 


0.009 


-1.135 


-1.100 


0.011 


-1.532 


-1.477 


0.013 


-1.444 


-1.420 


0.015 


-0.523 


-0.517 


0.017 


1.685 


1.695 


0.019 


5.667 


5.697 



TABLE VI: As Table |Vl but for spin polarized liquid ^He at ^ = 2/3. 



Density (A"^) 


Microscopic Calculation 


Parabolic Approximation 


0.003 


0.069 


0.124 


0.005 


-0.279 


-0.164 


0.007 


-0.652 


-0.488 


0.009 


-1.018 


-0.867 


0.011 


-1.382 


-1.164 


0.013 


-1.283 


-1.186 


0.015 


-0.349 


-0.331 


0.017 


1.874 


1.895 


0.019 


5.868 


5.891 
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FIG. 1: Two-body correlation function of liquid ^He, /(r), as a function of interatomic distance 
(r) at different values of spin asymmetry parameter ^ = 0.0, 1/3, 2/3 and 1.0 for p = 0.01 A~^. 
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FIG. 2: Internal energy of liquid ^He versus total number density at ^ = 0.0 (full curve), 1/3 
(dashed curve), 2/3 (dotted curve) and 1.0 (heavy dotted curve). 
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FIG. 3: Comparison of our results for the unpolarized (UnPol) and fully polarized (Pol) cases with 
the results of GFMC for different choices of the wave function. J refers to Jastrow, and J+T refers 
to Jastrow plus three-body wave functions [l^ . 
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FIG. 4: Comparison of our results for the unpolarized (UnPol) and fully polarized (Pol) cases with 
the results of FHNC [10]. J refers to Jastrow wave function. 
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FIG. 5: Comparison of our results for the unpolarized (UnPol) and fully polarized (Pol) cases 
with the results of FHNC [l^. J+T refers to Jastrow plus three-body wave function. 
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FIG. 6: Comparison of our results for the unpolarized (UnPol) and fully polarized (Pol) cases with 
the results of FHNC [10]. J+T+B refers to Jastrow plus three-body plus backflow wave function. 
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